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ABSTRACT 

The identification of a unique set of system parameters in large space struc- 
tures poses a significant new problem in control technology. This paper presents 
an infinite-dimensional identification scheme to determine system parameters in 
large flexible structures in space. The method retains the distributed nature of 
the structure throughout the development of the algorithm and a finite-element 
approximation is used only to implement the algorithm. This approach elimi- 
nates many problems associated with model truncation used in other methods of 
identification. The identification problem is formulated in Hilbert space and an 
optimal control technique is used to minimize weighted least squares of error be- 
tween the actual and the model data. A variational approach is used to solve the 
problem. A costate equation, gradients of parameter variations and conditions for 
optimal estimates are obtained. Computer simulation studies are conducted us- 
ing a shuttle-attached antenna configuration, more popularly known as the Space 
Control Laboratory Experiment (SCOLE) as an example. Numerical results show 
a close match between the estimated and true values of the parameters. 


PRECEDING PAGE BLANK NOT FILMED 


17 



DISTRIBUTED PARAMETER IDENTIFICATION 


TWO APPROACHES i 

FINITE-DIME SION AL METHOD 
INFINITE-DIMENSIONAL METHOD 



Table 1 


Literature Surveyed on the Parameter Estimation of Large Space Structures. 


Reference 

Approach 

Wells and Spalding (1977) [2] 

A finite-dimensional design approach 

Tung (1981) [3] 

where the structural model is trunc- 

Balas and Lilly (1981) [4] 

ated and the estimator is designed 

Balas (1981) [5] 

based on the reduced-order model. 

Lee and Bitter (1981) [6] 


Banks (1982) [7] 


Hendricks et al (1982) [8] 


Hendricks et al (1984) [9] 


Banks and Rosen (1984) [10] 


Rajaram and Junkins (1985) [ 1 1 J 


Lee, Walker and Hossain (1985) [12] 


Lee (1986) [13] 


Spalding (1976) [14] 

An infinite-dimensional design appr- 

Burns and Cliff (1977)[15[ 

oach where the PDE model is retained 

Sun and Juang (1982) [16] 

as long as possible and truncation is 

Lee (1986) [13] 

carried out only after the estimation 


algorithm is developed. 
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I World I Conceptual World 




I 


SYSTEM IDENTIFICATION SCHEME 




A Distributed System Model 


m(x)^(x,<) + D 0 -^{x,t) + A 0 u(x,t ) = F(x,t), 


xen, te{o,T\, 


F[x,t) = F B [x,t) + F c (xj) + F v (x,t), 


M 

Fc(x,0 = £o/ = X>(z)/.(0. 

t=l 


N 

F B {x,t) = B b g = ^6,(x)y,(i), 

i=i 


y = Cqu + ^out, 
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Basic Problem Formulation 


7 * 

A<) = ^ / o (y - z) T R{t){y - z)dt, 


where z is the measurement of output vector y given as 

z{x,t) = y(i, t) + e(x,f) 

with a measurement error e{x,t). Also, it is defined that 

{y - z) T R(t)(y - z) = f (y(x,f) - z{x,t)] T R{x,t){y{x,t) - z{x,t)\ dx, 

J n 
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Infinite-Dimensional Formulation 


£.<«) + + A(,M<) = «(»)/(<) »n(0, r|, 

u 6I 2 (0J;V), 6L 2 (0,T;/f), 

where /(f) is given in L 2 ( 0 ,T; V), and the initial conditions are 

u(0) = u 0 , u 0 given in V, 

and 

Q 

— u(0) = ui, U! given in H. 

dt 

The output function is 

y(0 = Cu (0» 

The identification problem can now be formulated as an abstract problem of deter- 
mining the parameter vector <7*(i) 6 Q that minimizes 

J(i ) = ^ /jy(0 - *(01 T *(0W0 - *MI*. 

where z(t ) is the observed data belonging to Y 
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Development of Infinite-Dimensional Identfication Algorithm 


THEOREM : Given a state equation (18) with initial conditions given by Eq. (19) and 
the cost function by Eq. (22) with y(t) satisfying Eq. (21), then the optimal parameter 
vector q satisfies the state equations (18)-(19) and the foilwing system of equations : 

£p(0 - D ■ jAt) + A- Pit) = -k T /!(Cu - *), (23) 

at- at * 

with the final conditions 

p(T) = j t v{T) = 0. (24) 

and the first variation of an augmented cost functional is 

SJ a = [ p T ^-{D^ + Au - Bf}6qdt = 0, (25) 

Jo <?<7 dt 

where p(t) is a costate variable also belonging to the Hilbert space V . 

PROOF : By combining Eqs. (18) and (22) an augmented cost functional can be defined 
as 

T 

J«(») = ^ f (»(0 - z(t)\ T Rit)\y(t) - z(t )\dt 

+ f T p(<f{£i “(0 + + Au[l)-Bf(t)] dt. (26) 
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Parameter Identification of Vibrating Beams 


Case I : A Simply-Supported Beam 

d * Q 4 

pA ~d^ + E1 d^ = X € [0,1), t > 0, 


“( x *0 = = °. x€d[0,L), < > o, 

tt(x,0) = - u (x,()| i=o = 0 , x€|0,£|. 


»(<) = 

2 (0 =y(0 + «(<)• 


j = 27 jf l» - *f% - *] *• 
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or 


d 2 u 

El d 4 u 1 

m 2 = 

m dx 4 m 

d 2 u 

d A u 

dt 2 ~ 



f) 


where 


9i 


EJ_ 

m 



and the parameter vector is defined by q — |gi,g 2 ] r • 


W- = f? + f 1“ ■ z|<(l " f >' 1 S [0, I|, i 6 [0,7% 

P( X > T ) = af p ( z ’*)lt-r = 0 ’ *G(0,I]. 

p(*,0 = ^ 2 P( I ’ f ) = °» 16 d(o,L],t e (o,r). 


_ f T ^pd 2 U 

6 gi “ y 0 5x2 5l 2 dt ' 

6J a f T L 

nr-L 

Thus, parameters can be updated by the steepest descent algorithm 
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Table 2 


Performance Data for Case I 


Iteration 

< 7i 

92 

m 

El 

/ 0 T error 2 dt 

1 

301.58 

0.0158 

63.00 

19000.00 

0.20178E-01 

2 

311.68 

0.0154 

64.84 

20207.72 

0.16640E-01 

3 

325.65 

0.0151 

66.27 

21582.87 

0.71528E-02 

4 

338.48 

0.0149 

66.91 

22647.68 ~ 

0.57208E-03 

5 

342.11 

0.0149 

66.98 

22913.64 

0.34507E-04 

6 

342.96 

' 0.0149 

66.99 

22973.13 

0.26524E-04 

7 

343.19 

0.0149 

66.99 

22989.35 

0.24364 E-04 

True 






values 

343.28 

0.0149 

67.00 

23000.00 
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Case II : A Cantileverd Beam 


m 





x € [0, L], t > 0, 


where £ is a damping coefficient. 


y(0 = u{L,t). 


The boundary conditions for a cantilevered beam are : 


d 2 

«(0,t) = 


= 0, 


i=0 


d 2 

dx 2 


u{x,t) 


z — L 


= Tz 3“^-0 


= 0 , 


z — L 


t > o, 

t > 0. 


The beam is initially at rest and hence the initial conditions are 


*(*.°) = = °> x6[0,I]. 


du(x, t) 


dx 


t = o 


dx 2 


u(x ,t) 


= 0, x e [o, L\. 


t = 0 


where 


d 2 u _ d z u d 4 u 

Tt r = ,3 5?ai " q ' aP + ' i)/( ' K 


El 

9 1 = — , 
m 

1 

?2 = — , 

m 

93 = 

V m 


and the parameter vector is defined by 9 = {91,92, to] 7 
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DISPLACEMENT 



TIME (ma) 
(b) 


Fig. 3 (a) Cantilevered beam with a step load, 
(b) resultant displacements. 
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Table 3 


Performance Data for Case II 


Iteration 

?i 

92 

93 

m 

El f 

Jq error 2 c it 

1 

301.59 

0.0159 

2.08 

63.00 

19000.00 0.060 

0.7843E-01 

2 

318.01 

0.0155 

2.76 

64.50 

20512.92 0.077 

0.1664E-01 

3 

326.20 

0.0153 

3.06 

65.26 

21286.34 0.085 

0.9432E-02 

4 

331.30 

0.0152 

3.23 

65.73 

21775.47 0.089 

0.4498 E-02 

5 

334.72 

0.0151 

3.34 

66.05 

22107.25 0.091 

0.2284 E-02 

6 

337.10 

0.0151 

3.43 

66.27 

22339.85 0.093 

0.1238E-02 

7 

338.83 

0.0151 

3.49 

66.43 

22510.09 0.095 

0.6649E-03 

8 

340.07 

0.0150 

3.53 

66.55 

22632.61 0.096 

0.3796E-03 

9 

340.96 

0.0150 

3.55 

66.65 

22721.78 0.096 

0.2199E-03 

10 

341.56 

0.0150 

3.56 

66.71 

22783.93 0.096 

0.1519E-03 

11 

342.06- 

0.0150 

3.59 

66.76 

22835.96 0.097 

0.8667E-04 

12 

342.44 

0.01497 

3.62 

66.80 

22874.19 0.098 

0.4961 E-04 

13 

342.73 

0.01496 

3.64 

66.83 

22903.95 0.098 

0.1165E-04 

True 

values 

343.28 

0.0149 

3.70 

67.00 

23000.00 0.1 
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Case III : A Simply-Supported Beam with Spatially Variable Parameter 


d 2 u 





b{x)f{t), 


where q = EI[x). 



d*_ 

dx 2 




^a[l) _ f T d^p d 2 U 

6q J 0 dx 2 dx 2 


x € [0, L], t > 0, 


i 6 [0,1], t€ [0,7), 
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SECTION 1 


STEP LOAD 


SECTION 2 


£ 


h L/2 

SECTION 1 
m = 67 k g / m 
EI = 23000 N/m 2 
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F>g- 4 (a) Simply-supported beam with spatially variable flexible rigidity, 

(b) resultant displacements at L/4, 

(c) resultant displcements at 3L/4. 
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Table 4 


Performance Data for Case III 


Iteration 

q{ sec. 1) 

<?(sec. 2) 

/ 0 T error 2 dt 

1 

21000.00 

39000.00 

0.76863E-02 

2 

21686.15 

39286.47 

0.40601E-02 

3 

22261.30 

39521.45 

0.14636E-02 

4 

22630.66 

39688.41 

0.41442E-03 

5 

22814.49 

39759.51 

0.13629E-03 

6 

22911.04 

39793.95 

0.51151E-04 

7 

22968.10 

39815.07 

0.1 1209 E-04 

True 




values 

23000.00 

40000.00 
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Fig. 1 The shuttle / antenna configuration of the spacecraft 
control laboratory experiment ( SCOLE ). 





Fig. 2 Drawing showing the direction of “roll" bending, 
“pitch" bending and “yaw" twisting. 
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The SCOLE Mathematical Model 


A. Dynamic Equations 


Roll Beam Bending Equation in v-z Plane 


+ EI *- d pr 


— f /^,n(0^(’ S 3 n) + 9<b,n{t) ^ > 

n=l ' 

t > 0, —00 < 3 < 00, 0 < 3 n < L. 


Pitch Beam Bending Equation in x-z Plane 


pA if ~ + EI, h ¥ 


— 3 n) + 98,n[t) g 3 ^( S ~ 3n )^ » 

t > 0, —CO < 3 < OO, 0<3 n <L. 


Yaw Beam Torsion Equation for z-Axis 


4 

~ ff^,n(t)<5(j — 3 n ), 

n=l 

t > 0, —00 < 3 < 00, 0 < 3 n < L. 
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B. Forcing Functions 


The forcing functions on the right side of each equation are dependent on boundary 
conditions and proof-mass actuators. 


Forces at 3 = si = 0 ( shuttle body forces ) 


The forces at s\ = 0 involves the shears at that point which are equal to the shuttle 
mass mi times the corresponding component of acceleration. 


UA*) = 

d 2 

= -m i — u*(0,t). 


Forces at 3 — 3 a — L ( reflector body forces ) 

d ^ 

d 2 d 2 

ft, 4(0 = -m A — u e (L,t) + m 4 r y — u+(L,t) + F x , 


where m 4 is the reflector mass, (r r ,r y ) is center of reflector mass from the beam tip at 
s = L, and F z and F v are the applied forces at the center of the reflector mass. 


Forces at s = s-> ( proof-mass actuator forces ) 


UAO = +m 2 — A^ i2 , 

ft, 2(0 = -m 2 — u e [s 2 ,t) + m 2 — A fl>2 , 


where A and m denote displacement and mass of the proof-mass actuator. 


Forces at 3 = 33 f proof-mass actuator force ) 



C. Moments 


Moments at s = 0 ( shuttle body moments ) 


(9+ A 
I 9e, 1 

\9i>,iJ 


— + wx ® IiWi] + Mi[t) + Mo{t), 


where I\ b the moment of inertia of the shuttle body, Afi(t) and M,o(t) are control and 
disturbance moments, respectively, applied to the shuttle body, and ® denotes the vector 
product. 

Moments at s = L f reflector body moments 1 


g *A \ /. . \ a 2 £4 

98,4 I = ~ M4««4 + «>4 ® Uu >4 ~ A/ 4 (0 + r ® F 4 (t)J - m 4 r ® 

9*. 4 J X ' 

where M 4 and F 4 are the control moment and force applied at the reflector center of the 

mass and £4 is the coordinates of the beam tip. 

Abo, J4 b the moment of inertia of the reflector, and I 4 b that with respect to the beam 
tip given by 

( r l - r » r » 0 \ 

I 4 = ^4 + ™4 ~ T z T y r z 0 

VO 0 rl + rjj 




Abstract Formulation of the SCOLE Problem 


M 0 r{t) + A 0 r{t) + B 0 F(t) + K 0 {r{t) 2 ) = 0 , 


where Mq is the 17 x 17 matrix specified by 


*i 

*2 

*3 

*4 

*5 

*0 

*T 



*4 
** 
*10 
*11 
*13 
*13 
*14 
*1$ 
*10 
*1T \ 


*3 

0 

pA 

0 


*3 

0 

0 

pA 


*4 


mi 

0 


*6 *4 


*7 *• *« *10 *11 



*13 *13 *14 *15 *10 *17 


0 m*r. 

0 -m«r f 
0 0 
0 0 
0 0 

/« 

m j 0 
0 mj 

m 3 0 
0 m 3 ; 
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Table 2 


Performance Data for Case I : Nonlinear SCOLE Model 


Iteration 

9 1 

92 

m 

El 

Jq error 2 dt 

1 

4.2353E+08 

11.765 

0.0850 

3.6000E+07 

0.1341E-04 

2 

4.2531E+08 

11.067 

0.0904 

3.8427E+07 

0.1896 E-05 

3 

4.2612E+08 

10.809 

0.0925 

3.9422E+07 

0.3038E-06 

4 

4.2743E+08 

10.705 

0.0934 

3.9929E+07 

0.9131E-07 

5 

4.2704E+08 

10.665 

0.0938 

4.0041E+07 

0.8031E-07 

6 

4.2617E+08 

10.637 

0.0940 

4.0064E+07 

0.6646E-07 

7 

4.2508E+08 

10.611 

0.0942 

4.0059E+07 

0.4938E-07 

8 

4.2390E+08 

10.589 

0.0944 

4.0029E+07 

0.3237E-07 

9 

4.2277E+08 

10.574 

0.0945 

3.9982E+07 

0.2032E-07 

10 

4.2174E+08 

10.561 

0.0947 

3.9934E+07 

0.1465E-07 

11 

4.2072E+08 

10.546 

0.0948 

3.9894E+07 

0.1332E-07 

12 

4.1959E+08 

10.525 

0.0950 

3.9867E+07 

0.1365 E-07 

13 

4.1819E+08 

10.495 

0.0953 

3.9845E+07 

0.1676E-07 

14 

4.1783E+08 

10.471 

0.0955 

3.9903E+07 

0.7113E-08 

True 






values 

4.1858E+08 

10.465 

0.0956 

4.0000E+07.00 
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Table 3 


Performance Data for Case I : Linearized SCOLE Model 


Iteration 



m 

El 

Jq error 2 dt 

1 

4.2353E+08 

11.765 

0.0850 

3.6000E+07 

0.1341E-04 

2 

4.2531E+08 

11.068 

0.0904 

3.8427E+07 

0.1894E-05 

3 

4.2613E+08 

10.809 

0.0925 

3.9422E+07 

0.3051E-06 

4 

4.2744E+08 

10.704 

0.0934 

3.9929E+07 

0.9142E-07 

5 

4.2705E+08 

10.665 

0.0938 

4.0041E+07 

0.8040E-07 

6 

4.2618E+08 

10.637 

0.0940 

4.0064E+07 

0.6660 E-07 

7 

4.2509E+08 

10.611 

0.0942 

4.0060E+07 

0.4955E-07 

8 

4.2392E+08 

10.590 

0.0944 

4.0030E+07 

0.3255E-07 

9 

4.2277E+08 

10.574 

0.0945 

3.9982E+07 

0.2032E-07 

10 

4.2175E+08 

10.561 

0.0947 

3.9933E+07 

0.1478E-07 

11 

4.2066E+08 

10.546 

0.0948 

3.9888E+07 

0.1378E-07 

12 

4.1957E+08 

10.524 

0.0950 

3.Q567E+07 

0.1361E-07 

13 

4.1816E+08 

10.495 

0.0953 

3.9846E+07 

0.1667E-07 

14 

4.1773E+08 

10.471 

0.0955 

3.9893E+07 

0.8371E-08 

True 






values 

4.1858E+08 

10.465 

0.0956 

4.0000E+07.00 
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FLEXURAL RIGIDITY 



0 2 4 6 8 10 12 1 4 

ITERATION 

(a) 



0 2 4 6 8 10 12 14 

ITERATION 

(b) 


Fig. 5 Convergence of parameters for the SCOLE problem in 
Case I using nonlinear model. 

(a) For pA / unit length. 

(b) For flexural rigidity, El. 
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>0 


TIME (■•) 


Fig- 4 Roll, pitch and yaw displacements when a damping 


lactor of 0.003 is added to the system. 







Table 4 


Performance Data for Case II 


Iteration 

El 

pA 

pl$ 

€ 

1 

32000000.00 

0.0750 

1.1000 

0.00250 

2 

37582208.00 

0.0881 

0.9690 

0.00250 

3 

39607652.00 

0.0928 

0.9428 

0.00250 

4 

40157244.00 

0.0941 

0.9362 

0.00250 

5 

40372424.00 

0.0946 

0.9247 

0.00280 

6 

40338068.00 

0.0945 

0.9183 

0.00288 

7 

40309128.00 

0.0945 

0.9127 

0.00293 

8 

40282296.00 

0.0943 

0.9088 

0.00289 

9 

40260868.00 

0.0943 

0.9047 

0.00287 

10 

40243016.00 

0.0942 

0.9008 

0.00285 

11 

40226244.00 

0.0942 

0.8983 

0.00281 

12 

40202376.00 

0.0942 

0.8979 

0.00286 

13 

40181576.00 

0.0942 

0.8977 

0.00292 

• 

14 

40162324.00 

0.0941 

0.8974 

0.00299 

15 

40149520.00 

0.0941 

0.8973 

0.00298 

16 

40140040.00 

0.0941 

0.8973 

0.00295 

True 





values 

40000000.00 

0.0956 

0.9089 

0.003 
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0 2 4 6 8101214 

ITERATION 

(a) 



0 24 6 8 10 12 14 

ITERATION 

(b) 


Fig. 7 Convergence of error criterion for Case II. 

t 

(a) For pA / unit length. 

(b) For pl^ /unit length. 
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CONCLUSION 


Infinite-dimensional identification method presented in this paper shows a significant 
promise in the parameter estimation of flexible structures with great potentials for appli- 
cations to LSS’s. The basic approach is the abstract formultion of the system dynamics 
in function spaces and then applying optimal control theory to adjust system parame- 
ters so that the error between actual and model data is minimized. The use of partial 
differential equation for the purpose of estimation elimainates many problems associated 
with model trunction in the finite dimensional approach. Based on partial differential 
equation models and a quadratic performance index an algorithm to estimate the opti- 
mal parameters has been developed. The numerical results show the effectiveness of the 
algorithm in estimating parameters of the flexible beam in the SCOLE problem. The 
results show fairly good match between the model and the estimated parameters. How- 
ever, as the number of parameters to be identified increases it becomes increasingly time 
consuming and difficult to solve. Also, due to model mismatch, slightly less accuracies 
are expected if experimental measurement data from physical beam were used. 
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